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Abstract: Time-series of marine inherent optical properties (IOPs) from 
ocean color satellite instruments provide valuable data records for studying 
long-term time changes in ocean ecosystems. Semi-analytical algorithms 
(SAAs) provide a common method for estimating IOPs from radiometric 
measurements of the marine light field. Most SAAs assign constant spectral 
values for seawater absorption and backscattering, assume spectral shape 
functions of the remaining constituent absorption and scattering 
components (e.g., phytoplankton, non-algal particles, and colored dissolved 
organic matter), and retrieve the magnitudes of each remaining constituent 
required to match the spectral distribution of measured radiances. Here, we 
explore the use of temperature- and salinity-dependent values for seawater 
backscattering in lieu of the constant spectrum currently employed by most 
SAAs. Our results suggest that use of temperature- and salinity-dependent 
seawater spectra elevate the SAA-derived particle backscattering, reduce 
the non-algal particles plus colored dissolved organic matter absorption, 
and leave the derived absorption by phytoplankton unchanged. 
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1. Introduction 

Satellite ocean color instruments provide consistent and high-density data records at sufficient 
temporal and spatial scales to allow retrospective analysis of long-term oceanographic trends. 
For example, the daily, synoptic images captured by the NASA Moderate Resolution Imaging 
Spectroradiometer onboard Aqua (MODISA) provide viable data records for observing 
decadal changes in biogeochemistry of both global and regional ecosystems [1]. Satellite 
ocean color instruments measure the spectral radiance emanating from the top of the 
atmosphere at discrete visible and infrared wavelengths. Atmospheric correction algorithms 
are applied to remove the contribution of the atmosphere from the total signal and produce 
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estimates of remote sensing reflectances (R rs (fy; sr _1 ), the light exiting the water normalized 
to a hypothetical condition of an overhead Sun and no atmosphere [2]. Bio-optical algorithms 
are applied to the R rs QC) to produce estimates of additional geophysical properties, such as 
spectral marine inherent optical properties (IOPs), namely the absorption and scattering 
properties of seawater and its particulate and dissolved constituents [3,4]. Time-series of these 
geophysical properties provide unparalleled resources for studying carbon stocks, 
phytoplankton population diversity and succession, and ecosystem responses to climatic 
disturbances on regional to global scales (e.g., [5-8]). 

Semi-analytical algorithms (SAAs) provide one mechanism for estimating marine IOPs 
from R rs (k), through a combination of empiricism and radiative transfer theory. Many 
published SAAs attempt to simultaneously estimate the magnitudes of spectral backscattering 
by particles (b bp (k); m -1 ), absorption by phytoplankton (a ph (X); m -1 ), and the combined 
absorption by non-algal particles and colored dissolved organic material (a dg (X); m -1 ) [3,4 and 
references therein]. These SAAs generally assign constant spectral values for seawater 
absorption and backscattering (a w (X) and b bw (X), respectively; m -1 ) [9,10], assume spectral 
shape functions of the remaining constituent absorption and scattering components, and 
retrieve the magnitudes of each remaining constituent required to match the spectral 
distribution of R r s(fy- As most SAAs differ only in the assumptions employed to define 
component spectral shapes and mathematical methods applied to calculate the magnitudes, 
the community typically focuses solely on refining these varied assumptions and methods [3]. 

While it is well known that a w (k) and b bw (X) in the visible spectrum vary with ocean 
temperature ( T ; °C) and salinity (S; g kg -1 ) [11-16], the impact of using constant values in 
lieu of temperature- and salinity-dependent values in an SAA has yet to be fully explored or 
quantified. Analyses conducted by Morel [11] revealed scattering increases of 30% in the 
spectral range 366-578 nm for S = 38.4 g kg -1 compared to those for pure water. This 
dependency of b bw (k) on S was most recently revisited by Zhang [16], whose theoretical 
model agrees to within 1% on average with Morel [12] at S = 38.4 g kg -1 . Despite this 
evolving understanding of how the spectral shape and magnitude of b bw (X) depend on T and S , 
the ocean color community continues to exclusively adopt the invariant Morel b bw (X) values 
[9,11,12] for use in SAAs and other bio-optical modeling activities [3,4], even though these 
values apply strictly to a singular temperature and salinity combination. The Zhang model 
suggests that b bw (k,T,S) can deviate from b bw (X,20,38A) by up to 25% at 443 and 547 nm, 
respectively, for typical temperature and salinity ranges found in natural waters (-2 < T < 
40°C and 0 < S < 40 g kg -1 ) (Fig. 1). Given that b b w( ^) — b bp (X) for some wuvelengths m many 
offshore and open ocean waters [15], one might expect that including temperature and 
salinity-dependent b bw (X) in an SAA (in place of a constant spectrum) will spatially and 
temporally alter the magnitudes of its derived products. 

We executed a series of analytical experiments to evaluate this expectation. In this paper, 
we present a preliminary quantification of the changes realized in SAA-derived products 
when using temperature- and salinity-dependent b bw (k) in lieu of a constant spectrum. We 
generated this quantification by applying the default SAA described in Werdell [3] (GIOP) to 
a wide dynamic range of in situ data and a global MODISA time-series using both b bw (k) 
from Morel [12] (which is the default for GIOP) and from Zhang [16]. For the latter, we 
calculated b bw (k,T,S) using the NOAA World Ocean Atlas sea surface salinity climatologies 
[17] and the NOAA optimally interpolated sea surface temperature climatologies [18] as 
inputs. While our use of climatologies may not have perfectly captured the potential 
differences realized using real-time data, they effectively enabled these proof-of-concept 
analyses and indicated that significant relative biases exist between b bp (X ), a dg QC), and a ph QC) 
derived from the SAA runs with and without temperature- and salinity-dependent b bw QC). We 
conclude with a brief discussion regarding the use of temperature- and salinity-dependent 
a w (k) in an SAA, however, we do not quantify its impact. 
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Fig. 1. Ratios of b bw (X,T,S) to b bw (K 20,38.4) for 443 nm (left panel) and 547 nm (right panel) 
for the temperature and salinity ranges -2 < T < 40°C and 0 < S < 40 g kg -1 . 

2. Methods 


2.1 A general semi-analytical algorithm 


Ocean color satellite instruments provide estimates of R rs (X), which can be converted to their 
subsurface values using the method presented in Lee [19]: 


r„a, 0 -) 


0.52 + 1.7 R rs (A)' 


( 1 ) 


Subsurface remote-sensing reflectances relate to marine IOPs via: 


r„(A, 0") = gj u(A) + g 2 u(A) 2 


u(A) = 


bM) 

a{X) + b b {X) 


( 2 ) 


where b b is the total backscattering coefficient (m _1 ), a is the total absorption coefficient 
(m _1 ), and g (sr -1 ) vary with illumination conditions, sea surface properties, and the shape of 
the marine volume scattering function. GIOP defaults to g 7 = 0.0949 and g 2 = 0.0794 sr -1 
from Gordon [20]. The absorption coefficient can be expanded as the sum of all absorbing 
components. Further, each component can be expressed as the product of its concentration- 
specific absorption spectrum (eigenvector; a*) and its magnitude (eigenvalue, 47): 

a(l) = a w (l) + M ph a ph (l) + M dg a dg (l\ (3) 

GIOP expresses a* dg (k ) as exp(- s dg X), where $dg describes a rate of exponential decay (nm l ), 
and assigns $dg 0 .018 nm '. The term M ph is equivalent to the concentration of the 
phytoplankton pigment chlorophyll-a ( C a ; mg m -3 ) when a * p h(X) is expressed as a 
chlorophyll-a specific absorption spectrum (m 2 mg -1 ). GIOP uses the chlorophyll-specific 
a* ph (k) from Bricaud [21], seeded using an estimate of C a from O’Reilly [22] and normalized 
to a* p h( 443) = 0.055 m 2 mg -1 . Similar to a(X ), b b (X) can be expanded to: 

b b (l) - b bw (l) + M bp b* bp (l). (4) 

GIOP expresses b* bp (X) as X Sbp , where S bp defines the steepness of the power law (we 
acknowledge the validity of this power function remains debatable), and estimates S bp using 
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the R rs (X) band-ratio approach described in Lee [19]. Werdell [3] describes in detail the 
rationale for the selection of the eigenvectors used in Eqs. (3) and (4). 

Using R rs (k ) and eigenvectors as input, the eigenvalues M can be estimated via linear or 
nonlinear least squares inversion of Eqs. (1) to (4). GIOP employs the nonlinear Levenberg- 
Marquardt method and uses all R rs (X ) between 400 and 700 nm. We considered retrievals of 
b bp (X), a dg (k), and a ph (k) to be valid when they fell within predefined ranges (e.g., -0.05 a w (X) 
< a dg (X) < 5 nT 1 ) and could reconstruct R rs (X ) that differed from the input spectrum by less 
than 33% anywhere in the range 400-600 nm. Werdell [3] describes the additional algorithm 
metrics and quality control practices embedded within GIOP. 

2.2 Data acquisition 

We acquired coincident observations of in situ R rs (X), b b (X), a dg (X), and a ph (X) from the NASA 
bio-Optical Marine Algorithm Data set (NOMAD) [23]. Sample sizes for b b (X) ( N = 217) 
differ from a dg (X) and a ph (k) (N= 609) in NOMAD because of instrumental variability in field 
campaigns included in its compilation. We obtained MODISA Level-3 R rs (X) monthly 
composites for 2010 from the NASA Ocean Biology Processing Group (OBPG; 
http://oceancolor.gsfc.nasa.gov). We also acquired NOAA World Ocean Atlas 2009 sea 
surface salinity monthly climatologies (WOA-SSS) [17] and NOAA optimally interpolated 
sea surface temperature weekly climatologies (OISST) [18] from the OBPG, as they are 
available through the NASA SeaWiFS Data Analysis System (SeaDAS; 
http://seadas.gsfc.nasa.gov). The NOAA National Oceanographic Data Center and NOAA 
Earth System Research Laboratory provide and maintain the official source WOA-SSS and 
OISST data records, respectively. 

Table 1. Ordinary linear regression statistics for modeled (GIOP) versus in situ 
(NOMAD) IOPs at MODISA visible wavelengths using b bw c (K) and b bw TS (k ) 


GIOP-C GIOP-TS 




r 2 

Slope ( ± SE) 

Bias* 

MPD + 

r 2 

Slope ( ± SE) 

Bias* 

MPD + 


412 

0.56 

0.88 (0.04) 

0.128 

27.8 

0.56 

0.98 (0.05) 

0.075 

20.6 


443 

0.60 

0.90 (0.04) 

0.132 

28.5 

0.60 

0.97 (0.04) 

0.091 

22.6 

bbp 

488 

0.65 

0.94 (0.04) 

0.121 

29.9 

0.65 

0.98 (0.04) 

0.098 

26.3 

531 

0.66 

0.95 (0.04) 

0.123 

31.7 

0.66 

0.97 (0.04) 

0.108 

27.7 


547 

0.69 

0.97 (0.04) 

0.119 

32.3 

0.69 

0.97 (0.04) 

0.110 

30.1 


667 

0.73 

1.00 (0.04) 

0.107 

30.9 

0.73 

0.98 (0.04) 

0.112 

31.2 


412 

0.76 

1.11 (0.02) 

-0.110 

28.6 

0.74 

1.13(0.03) 

-0.138 

31.6 


443 

0.72 

1.08 (0.03) 

-0.142 

35.0 

0.70 

1.10(0.03) 

-0.170 

36.6 

a dg 

488 

0.66 

1.03 (0.03) 

-0.212 

44.3 

0.64 

1.05 (0.03) 

-0.240 

45.3 


510 

0.64 

1.01 (0.03) 

-0.235 

46.7 

0.62 

1.03 (0.03) 

-0.264 

47.9 


547 

0.59 

0.98 (0.03) 

-0.303 

53.1 

0.57 

0.99 (0.03) 

-0.331 

55.3 


667 

0.46 

0.90 (0.02) 

-0.499 

69.4 

0.44 

0.92 (0.02) 

-0.526 

70.7 


412 

0.65 

1.14(0.03) 

-0.053 

25.3 

0.72 

1.09 (0.02) 

-0.050 

25.5 


443 

0.65 

1.14(0.03) 

-0.033 

23.2 

0.73 

1.08 (0.02) 

-0.030 

23.2 

Clph 

488 

0.66 

1.13(0.03) 

-0.003 

24.9 

0.74 

1.08 (0.02) 

0.000 

24.2 


531 

0.70 

1.14(0.03) 

0.008 

27.7 

0.77 

1.10(0.02) 

0.011 

26.8 


547 

0.69 

1.08 (0.02) 

0.073 

38.1 

0.74 

1.05 (0.02) 

0.077 

37.2 


667 

0.75 

1.04 (0.02) 

0.107 

34.8 

0.80 

1.01 (0.02) 

0.111 

34.4 


* Bias = X(model - in situ) / N, where N is the sample size. 
+ MPD = median(100% * |model / in situ - 1 1). 


2.3 Backs cattering of seawater 

As suggested above, most SAAs employ the b bw (X) provided by Morel [12], as tabulated via a 
best-fit expression by Smith and Baker [9]. This fit is commonly expressed as a power-law, 
namely b bw (k) = 0.0038 (400 / X)' 432 , where ft /m (400) = 0.0038 m~‘ (e.g., [24]). This 
expression provides the constant b bw (k ) used in this activity, which we hereafter refer to as 
b bw c (X). Using the WOA-SSS and OISST data records as input into Zhang [16] (via 
MathWorks MATLAB software kindly provided by X. Zhang), we calculated temperature - 
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and salinity-dependent b bw (X) for each station in NOMAD and pixel in the MODISA monthly 
composites. We hereafter refer to temperature- and salinity-dependent b bw Q2) as b bw TS (X). 

2.4 Data processing and analysis 

We applied GIOP to each station in NOMAD and each pixel in the MODISA monthly 
composites twice - first using b bw c (X) and second using b bw TS (X ), with all other parameters 
held constant. We hereafter refer to these two instances of GIOP as GIOP-C and GIOP-TS, 
respectively. We focused on two biogeochemically distinct regions in the MODISA monthly 
composites - the North Atlantic (40 to 50°N, -50 to -20°W) and South Pacific Gyre (-20 to 
-30°S, 100 to 130°E). We used OBPG-supported satellite data processing software (13 gen; 
distributed with SeaDAS) to run GIOP on the Level-3 MODISA R rs (X). Currently, 13gen 
accepts WOA-SSS and OISST data records as standard ancillary input. For NOMAD, we 
compared modeled and in situ b bp (k ), a dg (X), and a ph (X) for the two GIOP runs using ordinary 
least squares regression. Our statistics of interest included the coefficient of determination 
(r 2 ), the regression slope, the bias, and the absolute median percent difference (MPD; see 
Table 1 for its calculation). We executed these comparisons using in situ b bp (k) calculated 
with the corresponding b bw (X) from the GIOP run (i.e., b bp (k) = b b (X) - b bw c (X) for GIOP-C 
and b bp (X) = b b (k) - b bw TS (k ) for GIOP-TS). For both NOMAD and MODISA, we calculated 
the ratios of IOP retrievals from the two runs as GIOP-TS over GIOP-C. 








in situ b bp (443) 


in situ ^(443) 


in situ £^(443) 


Fig. 2. Modeled (GIOP) versus in situ (NOMAD) IOPs using b bw c (X) (A-C) and b bw TS (k) (D-F). 

The solid line shows a 1:1 relationship. 

3. Results 

Direct comparisons of modeled and in situ b bp (k ), a dg (k ), and a ph (X) provided estimates of the 
accuracy of GIOP-C and GIOP-TS. Visually, the model-versus -/>2 situ regression analyses for 
the two runs showed only subtle differences (Fig. 2). For MODISA wavelengths, GIOP-TS 
realized r 2 that were unchanged for b bp (X ), slightly degraded for a dg (X) (—3%), and modestly 
improved for a ph (k ) ( — F 10%) relative to GIOP-C (Table 1). GIOP-C showed spectral 
dependence in its regression slopes for b bp (k) (rising from 0.88 to 1.00 for 412 to 667 nm), but 
GIOP-TS did not, hovering around unity (0.97-0.98) for all wavelengths. GIOP-C and GIOP- 
TS yielded similar spectral dependence in the regression slopes for a dg (X) (falling from 1.11 to 
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0.90 and 1.13 to 0.92 for 412 to 667 nm, respectively). But, the regression slopes for a ph QC) 
fell closer to unity for GIOP-TS (average of 1.07) than for GIOP-C (average of 1.1 1). Finally, 
GIOP-TS realized improved (reduced) biases and MPDs for b hp (k) relative to GIOP-C, with 
the exception of b bp (661). GIOP-TS yielded slightly degraded biases and MPDs for a dg (k ,), but 
comparable statistics for a pb (K), relative to GIOP-C. 





0.90 0.95 1.00 1.05 1.10 0.90 0.95 1.00 1.05 1.10 0.90 0.95 1.00 1.05 1.10 

GIOP-TS / GIOP-C GIOP-TS / GIOP-C GIOP-TS / GIOP-C 

Fig. 3. Frequency of ratios of IOPs derived from GIOP-TS and GIOP-C run on NOMAD 
R rs QC)- The solid and dotted lines show unity and the distribution median. Bin sizes are 0.01. 

On average, estimates of b bp (k ) from GIOP-TS rose 2.7% compared to those from GIOP- 
C for the full population of data in NOMAD (Fig. 3). In contrast, estimates of a dg (k ) and 
a ph (k) fell 6.1 and 1.0%, respectively. These departures showed dependence on the trophic 
state of the station (e.g., low C a , oligotrophic water versus high C a , eutrophic water). In 
particular, GIOP-TS yielded much higher b hp (X) than GIOP-C (3-10%) for the clearest 
NOMAD stations ( C a < 1 mg m -3 ), whereas a dg (X) from GIOP-TS fell by 6% in this range 
(Fig. 4). Both GIOP-TS and GIOP-C produced comparable a ph (X ) over the full dynamic range 
of NOMAD C a , in particular when C a > 0.2 mg m -3 . GIOP-TS yielded slightly lower a ph { 443) 
than GIOP-C in the clearest waters, for example, where C a ~ 0.1 mg m -3 . 

MODISA global imagery from May 2010 showed geographic features that mimicked this 
“trophic-level” behavior (Fig. 5). In the oligotrophic sub-tropical gyres, for example, b bp { 443) 
from GIOP-TS exceeded that from GIOP-C by 10-20%, whereas a dg { 443) and a ph ( 443) fell 5- 
10% and 1-5%, respectively. By far, b bp ( 443) realized the largest spatial differences in 
magnitude, as most open ocean waters showed significantly elevated GIOP-TS to GIOP-C 
ratios. The Southern Ocean (circumpolar Antarctica) and near-shore waters presented the only 
exceptions, where b bp ( 443) was largely equivalent for both GIOP runs. The modest 
differences in a ph ( 443) followed a similar and equally broad spatial pattern to b bp ( 443), 
differing by a few percent in the open ocean and approaching unity in the Southern Ocean and 
along the coasts. In contrast, a dg { 443) from GIOP-TS and GIOP-C only truly differed in the 
centers of the sub-tropical gyres (and occasionally along the coast). While the Southern 
Ocean and coasts often realize the largest departures in T and S from the 20°C and 38.4 g kg -1 
inherent to b bw c (k) (Fig. 6), the differences between b bw TS (X) and b bw c (X) remain minimal for 
the T and S combinations of these regions (Fig. 1). 

Time-series of MODISA monthly averages in the North Atlantic showed seasonal patterns 
in the ratios of GIOP-TS to GIOP-C 6^(443) (Fig. 7). GIOP-TS produced b bp ( 443) that 
exceeded that from GIOP-C by 3-5% in April- July and by 7-12% in January-March and 
August-December. Ratios of GIOP-TS to GIOP-C 0^(443) followed a subtle, inverse 
seasonal pattern, although the two retrievals never departed by more than 3% (in boreal 
winter). The b bp ( 443) and a^(443) time-series tracked each other through a natural, seasonal 
progression of phytoplankton biomass [25], with higher b bp ( 443) generally occurring in the 
presence of higher a ph ( 443) (e.g., during and following the spring bloom). The highest GIOP- 
TS to GIOP-C ratios for b bp ( 443) appeared with the lowest b bp ( 443), while the highest ratios 
for a pb ( 443) appeared with the highest a ph ( 443) in the late boreal spring. Interestingly, despite 
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a spike in magnitude in boreal spring and a dip in boreal winter, the ratios of a</ g (443) showed 
no temporal dependence, but rather remained fixed with GIOP-C exceeding GIOP-TS by 6%. 

MODISA monthly averages in the South Pacific Gyre showed little seasonal pattern in the 
ratios of GIOP-TS to GIOP-C (Fig. 7). This region encompasses the clearest known ocean 
water [15] and the annual patterns in IOP magnitudes followed expectations, with minor 
variability in b bp ( 443), a peak in a ph ( 443) in June and July (austral winter), and a subsequent 
rise in a dg ( 443) [26,27]. Overall, GIOP-TS produced b bp { 443) that exceeded GIOP-C by 11%, 
a dg { 443) that fell below by 6%, and a ph ( 443) that fell below by 4% consistently over the full 
annual cycle. In general, the average magnitudes of these departures follow those observed 
for the clearest stations in NOMAD (e.g., C a < 0.1 mg m -3 ; Fig. 4). 
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Fig. 4. Ratios of IOPs from GIOP-TS and GIOP-C run on NOMAD R rs QC) as a function of their 
corresponding in situ C a . We considered ten equally spaced (logio) bins for the range 0.05 < C a 
<30 mg m“ 3 . The black circles show the median and are placed at the bin center. The solid 
vertical lines show the standard deviations for each C a bin. 

4. Discussion 

We demonstrated that using b bw TS (X ) (e.g., [16]) in an SAA in lieu of the community standard 
b b w C (fy [12] impacts the derived products in a sufficiently significant way as to merit further 
consideration. Theoretically and pragmatically, our results suggest that an SAA configuration 
such as GIOP-TS provides a step towards more accurate estimation of marine IOPs. 
Quantitatively, we expect our results could be refined with more contemporaneous values of 
T and S - for example, using coincident sea surface temperature data records from MODISA 
and complimentary sea surface salinity records from the Aquarius instrument onboard the 
Argentinean SAC-D spacecraft. Our case study cumulatively suggests that using b b J s (k) 
results in elevated b bp (k) (3-10%), reduced a dg (X) (1-6%), and negligibly changed a ph (k). 

Given that b bw TS (k) typically falls below b bw c (X) for most S<38.4 g kg -1 (Fig. 1), we 
expected GIOP-TS estimates of b bp (k) to exceed those from GIOP-C (Figs. 3-6). Also as 
expected, the largest departures in b bp (k) occurred in the open ocean where b bw ( X) > b bp (k) 
[15] (in contrast to highly turbid or productive environments where b bw (X) « b bp (k)). Of 
interest, however, which we elaborate upon later, is that b b (X) from GIOP-TS showed a 
consistently low bias relative to that from GIOP-C (Fig. 8). The NOMAD regression results 
indicate that b bp (k) from GIOP-TS realized improved agreement with in situ values (Fig. 2, 
Table 1). For the ocean color remote sensing paradigm, this implies that SAAs currently 
underestimate b bp (K ). Following, the switch to b bp TS (k ) could refine our future scientific 
understanding of marine particle dynamics and primary productivity (e.g., [7]). More 
compellingly, the MODISA time-series indicate that using b b J s (K) in place of b bw c (X) imparts 
an amplified seasonal signal in b bp (k ), which could further alter our understanding of and 
expectations for these phenomena. Finally, the spectral slope for b bw TS (X) also varies with 
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Fig. 5. MODISA Level-3 monthly IOPs for May 2010. Panels (A) to (C) show a p h( 443) from 
GIOP-C, a ph { 443) from GIOP-TS, and the ratio of the two (GIOP-TS / GIOP-C), respectively. 
Panels (D) to (F) show a^(443) from GIOP-C, <2^(443) from GIOP-TS, and the ratio of the 
two (GIOP-TS / GIOP-C), respectively. Panels (G) to (I) show b bp ( 443) from GIOP-C, 
bb P { 443) from GIOP-TS, and the ratio of the two (GIOP-TS / GIOP-C), respectively. 
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salinity (Zhang [16] reported a range of -4.286 to -4.306 from 0 to 40 g kg -1 ), whereas the 
slope of b bw c (X) does not, which will impact SAAs such as Loisel [5] that produce estimates of 
the spectral slope of b bp QC), S bp , to infer marine particle sizes. 


T - 20 (°C) S- 38.4 (a kg 1 ) 



-10 -8 -6 -4 -2 0 S-38.4 (cj, kg* 1 ) 


Fig. 6. OISST and WOA-SSS long-term climatologies for May, less 20°C (left) and 38.4 g 
kg -1 (right), respectively. 

Using b bw TS (k ) yielded reduced estimates of both a dg (k) and a ph (X ), although departures in 
the latter were largely negligible (Figs. 3-6). The NOMAD regression results indicate that 
a pb (X) from GIOP-TS improves upon that from GIOP-C, but that a dg (X) degrades slightly (Fig. 
2, Table 1). Werdell [3] demonstrated that SAA retrievals are highly sensitive to the 
parameterization of the spectral slope of a dg (X), S dg , and that changes in a dg (X) often 
accompany changes in b bp (X) (a classic case of transfer of variance, as a dg (X) and b bp (X ) have 
similar spectral shapes). However, repeating our experiment with S dg = 0.015 and = 0.02 
nm" 1 , and with Sdg dynamically estimated using the R rs (fy band-ratio approach described in 
Lee [19], produced equivalent results (not shown). Globally, the most significant departures 
occurred in subtropical gyres (Fig. 5). While this maintains some significance for studies 
relating to carbon budgets (particularly, the optically-relevant dissolved component, e.g., 
[6,8]), the bias between a dg (X) from GIOP-TS and GIOP-C remained fairly constant over most 
trophic levels (Fig. 4) and seasons (Fig. 7). Following, the switch to b bw TS (X) may produce 
a dg (k) with reduced magnitudes, but not with significantly altered spatial distributions. We 
recommend that future studies continue to focus on methods to dynamically assign S dg to each 
station in lieu of using spatially and temporally fixed values, with the purpose of reducing 
differences in SAA-derived and in situ a dg (K). 

Intuitively, the elevated estimates of b bp (k ) realized with b bw TS (k) appear in contradiction 
with the reduced estimates of a dg (k) and a ph (X). Given fixed a w (X) [10], these reduced 
estimates yielded a(X) from GIOP-TS that fell below those from GIOP-C by 3.2% on average 
(Fig. 8). The contradiction lies within Eq. (2) - that is, we expect common u(k) from both 
GIOP runs, assuming equivalent least-squares fitting performances, which suggests that b b (X) 
from GIOP-TS should be reduced in conjunction with a(X). With b bw TS (X) typically less than 
b bw c (fy, and b bp (k ) from GIOP-TS holistically greater than that from GIOP-C, we initially 
anticipated b b (k) from GIOP-TS to be roughly equivalent to that from GIOP-C. In other 
words, we expected a corresponding rise in magnitude in derived b bp (k) that accounted for the 
decline in magnitude in b bw TS (k). In practice, however, estimates of b b (X) from GIOP-TS fell 
3.5% below those from GIOP-C (Fig. 8). This decline in b b (X) ultimately yielded equivalent 
u(X) for the two runs (0.15% different on average). In the end, using b bw TS (k ) in lieu of b bw c (k ) 
yielded reduced b b (k ) and a(k ), despite producing elevated b bp (X). 

Changes in the parameterization of an SAA (e.g., the component spectral shapes) change 
the derived products to varying degrees [3]. Our results do not imply that GIOP-TS provides 
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the optimal SAA configuration for all water masses at all times [28]. As discussed in Werdell 
[3], significant work remains to refine SAAs for both global and regional applications. 
Westberry [29], for example, recently demonstrated that accounting for Raman effects can 
additionally improve SAA-estimates of b hp (k). Despite the work that remains to be 
accomplished, however, we felt it prudent to repeat several of our experiments with an 
alternative form of SAA, namely the Quasi- Analytical Algorithm (QAA [19];). A description 
of the differences in algorithm form relative to GIOP exceed the scope of this paper, but can 
be found in Werdell [3]. For the NOMAD R rs (X), using b bw TS (k ) in lieu of b bw c (X ) as input into 
QAA yielded elevated b bp (X) (3.2%), reduced a dg (k ) (5.5%), equivalent a ph (X ) (0.5%), reduced 
b b (X) (2.8%), and reduced a(X) (2.9%) on average. These departures effectively match those 
from GIOP in both magnitude and direction. 
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Fig. 7. MODISA Level-3 monthly time-series for 2010 in the North Atlantic and South Pacific 
Gyre. Black and red circles indicate GIOP-C and GIOP-TS, respectively. Blue circles how the 
ratio of the two (GIOP-TS / GIOP-C). 

Cumulatively, with the assumption that global climatologies of T and S provided 
sufficiently reliable values for this case study, our results imply that using temperature and 
salinity-dependent b bw (X) [16] in place of the (most commonly adopted) constant spectrum 
[12] in an SAA results in elevated b bp (X ) and reduced a dg (X). Spatially and temporally, a ph (k) 
from GIOP-TS and GIOP-C showed little variability, with the exception of the highest 
latitudes. We acknowledge that the departures in b bp (X) and a dg (X) may be sufficiently small to 
also be considered insignificant, however, it appears that using b bw TS (k ) produces improved 
comparisons with field data, at least for b bp (X) and a ph (k). Using b bw TS (X) in an SAA also 
appears to change the spatial and temporal distributions of b bp (X ), which should be 
sufficiently compelling itself to stimulate subsequent investigations, despite the degraded 
comparisons between modeled and in situ a dg (k). 
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Fig. 8. Frequency of ratios of IOPs and u( 443) from Eq. (2) derived from GIOP-TS and GIOP- 
C run on NOMAD R rs ( X). The solid and dotted lines show unity and the distribution median, 
respectively. Bin sizes are 0.01. 

We limited our activity to the use of b bw TS (X) in an SAA. In practice, b bw (X) appears 
elsewhere within the satellite ocean color paradigm, such as within the atmospheric correction 
process [2]. For example, the standard NASA algorithm applied to account for non-negligible 
near-infrared (NIR) radiances uses b bw (X) and a w (X) to model ^(NIR) [30]. Although this 
specific algorithm employs an iterative process, using b b J s (k ) in lieu of b bw c (k ) has potential 
to alter R rs ( NIP) and, thus, the aerosol models selected with the atmospheric correction 
process. Likewise, a w (X) appears throughout the ocean color paradigm. Unlike b bw (X ), the 
temperature- and salinity-dependence of a w (k) most strongly effects the NIR portion of the 
spectrum [13,14]. We, therefore, do not anticipate that using a w TS (X) in place of a w c (X) will 
significantly alter SAA-derived IOPs. With regards to the atmospheric correction component 
presented above, however, this replacement will almost certainly alter the modeled ^(NIP) 
derived following the method Bailey [30]. Accounting for non-negligible ^(NIP) remains a 
critical step within the atmospheric correction process in turbid waters ([31] and references 
therein). We recommend that similar, subsequent studies emerge to evaluate the use of 
temperature- and salinity-dependent a w (X) in standard satellite ocean color data processing. 
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